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We performed a systematic study on the properties of CeFeAsi-^P^O (0 < x < 1) by electrical 
resistivity, magnetization and specific heat measurements. The c-axis lattice constant decreases 
significantly with increasing P content, suggesting a remarkable chemical pressure. The Fe-3d elec- 
trons show the enhanced metallic behavior upon P-doping and undergo a magnetic quantum phase 
transition around x « 0.4. Meanwhile, the Ce-4/ electrons develop a ferromagnetic order near the 
same doping level. The ferromagnetic order is vanishingly small around x — 0.9. The data suggest a 
heavy-fermion-like behavior as x > 0.95. No superconductivity is observed down to 2 K. Our results 
show the ferromagnetic ordered state as an intermediate phase intruding between the antiferromag- 
netic bad metal and the nonmagnetic heavy fermion metal and support the cerium-containing iron 
pnictides as a unique layered Kondo lattice system. 

PACS numbers: 74.62.Dh, 74.62.Bf, 74.62.Fj, 71.20.Eh, 74.70.Dd 



I. INTRODUCTION 

The homologous quaternary iron pnictides, i?FeA"0, 
(R: rare earths, X: As or P) show a diversity of physical 
properties although they have the same ZrCuSiAs-type 
crystal structure with a unique FcA-laycr sandwiched by 
the #0-layeri. Typically, LaFeAsO is an antifcrromag- 
netic (AFM) metal below 140 K, and is a superconductor 
(SC) with maximal T c = 26 K upon F-doping^. T c can 
be surprisingly increased up to 41 K when La is replaced 
by Ce 3 - or even higher if replaced by other rare earths^—. 
Meanwhile, the parent compounds (denoted as Ln-1111) 
of these SC's, such as CcFcAsO, arc still Fe-3c? itiner- 
ant AFM's similar to LaFeAsO. However, a noticeable 
/-electron AFM ordering of Ce 3+ is also observed at a 
much lower temperature (~4 K) 3 -. On the other hand, 
the stoichiometric LaFePO is a low T c SC without any 
trace of the AFM ordering^. By contrast, CeFePO is 
a heavy fermion (HF) metal with Kondo temperature 
T K ~ 10 K2. 

These discoveries put the rare-earth iron pnictides on 
the boundary between the high-T c superconductors^ and 
the heavy fermion metals^, and open an avenue for 
searching the complicated interplay of various d- and /- 
electron correlations. Within the iron-pnictogen layer, 
the ci-electrons are expected to be more itinerant in the 
phosphides than in the arsenide s 12 ' 13 , with the former be- 
ing recently identified as a moderately correlated metalii. 
In contrast to the cuprates, the inter-layer distance be- 
tween the transition metals and the rare-earths in the 
iron pnictides is critically significant so that the cou- 
pling between them may play an important role even in 
the parent compounds of the high-T c Fe-based SC's^. 
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Recent neutron scattering and muon spin relaxation ex- 
periments indeed provide evidence for a sizable inter- 
layer coupling in CeFeAs0 16 i 17 . So far it is still unclear 
whether this coupling is due to some kind of polarization 
effect raised by the ordered moment of Fe 2+ -ions^ or 
more microscopically the effective hybridization between 
the 3d- and 4/-orbitals bridged by the pnictogensi^. A 
first-principle local-density approximation (LDA) plus 
dynamical mean- field theory (DMFT) studyiS. suggested 
that applying physical pressure will enhance the d-f hy- 
bridization, leading to the Kondo screening of the Ce- 
moments. However a sufficiently high pressure is re- 
quired in order to observe this effect. It turns out that 
CeFeAsi-^P^O, i.e., P doping at As sites in CeFeAsO, 
may provide, among others, a unique layered Kondo lat- 
tice system to probe the intriguing 3<i-4/ electron in- 
terplay under ambient pressure^. Chcmical-prcssurc- 
induced superconductivity has been observed in P-dopcd 
LaFeAsO^. 



We report a systematic study on the doping evolu- 
tion of the physical properties of CeFeAsi-^P^O using 
electrical resistivity p(T), magnetic susceptibility x(T), 
isothermal magnetization M(H) and specific heat C(T) 
measurements. A series of 21 different P-dopcd polycrys- 
tallinc samples were synthesized. Our results reveal a 
rich phase diagram consisting of an AFM quantum crit- 
ical point (QCP) of the Fe-3d electrons and a possible 
ferromagnetic (FM) QCP of the Ce-4/ electrons. An in- 
termediate FM phase emerges in between the AFM bad 
metal and the nonmagnetic HF metal. Because P dop- 
ing does not introduce extra electrons but shortens the 
c-axis significantly, these findings provide a rare example 
of chemical-prcssurc-induced HF metals with strong FM 
fluctuations. 



2 



II. EXPERIMENTAL 

We synthesized a series of CeFeAsi-^P^O (0< X < 
1) polycrystallinc samples by solid state reaction, during 
which, Ce, Fe, As, P and Ce0 2 of high purity(> 99.95%) 
were used as starting materials. First, CeAs (or CeP) 
was presynthesized by reacting Ce discs and As (or P) 
powders at 1323 K for 72 h. FeAs (or FeP) was prepared 
by reacting Fe and As (or P) powders at 1173 K (or 1023 
K) for 20 h. Second, powders of CeAs, CeP, Ce0 2 , FeAs 
and FeP were weighed according to the stoichiometric ra- 
tio, thoroughly ground, and pressed into a pellet under 
a pressure of 600 MPa in an Argon filled glove box. The 
pellet was sealed into an evacuated quartz tube, which 
was then slowly heated to 1448 K and kept at that tem- 
perature for 50 h. 

Powder X-ray diffraction (XRD) was performed at 
room temperature using a D/Max-rA diffractometer with 
Cu-K Q radiation and a graphite monochromator. Lat- 
tice parameters were refined by a least-squares fit us- 
ing at least 30 XRD peaks, and the structural refine- 
ments were performed using the programme RIETAN 
2000. The dc magnetization measurement was carried 
out in a Quantum Design Magnetic Property Measure- 
ment System (MPMS-5). Physical Property Measure- 
ment System(PPMS-9) was used to take the resistivity 
measurement, as well as specific heat. 



III. RESULTS AND DISCUSSION 

The room-temperature powder x-ray diffcraction 
(XRD) patterns of CeFeAsi„ x P :E are shown in Fig.l. 
All the XRD peaks can be well indexed based on the 
tetragonal ZrCuSiAs-type structure with the space group 
P4/nmm (No. 129), and no obvious impurity phases can 
be detected, suggesting high quality of the samples. As 
shown in Fig.l(a-b), the XRD peaks shift toward right- 
hand, especially for the high angle reflections. The (003) 
peaks shift much faster than the (110) peaks so that they 
merge into one peak in the high doping case, indicating 
that the c-axis shrinks more severely than the a-axis. 
This observation is consistent with the lattice parameter 
calculation results shown in Fig. 1(c). Fig.l(d-e) present 
the Rietveld refinement profiles for the two end members 
CeFeAsO and CeFePO, which give details of the crystal 
structure. The remarkable differences in structure lie in 
the positions of Ce and As/P. For CeFePO, the thickness 
of FeP layers is much smaller, meanwhile, Ce atomic lay- 
ers are much closer to the FeP layers. This structural 
feature supports the strong coupling between Fe-3d and 
Ce-4/ electrons. 

The temperature dependent resistivity is shown in 
Fig. 2, where four prominent features can be identified, (i) 
The overall suppression of the room temperature resistiv- 
ity with increasing doping. The resistivity at T — 300 K 
decreases gradually from 300 /iQm for x=0 to 12 [iQra for 
x=l, indicating that P-doping enhances the metallic be- 



havior significantly, (ii) The resistivity anomaly, mostly 
pronounced for x = 0. This anomaly was ascribed to 
the structure distortion and the accompanied Fe-AFM 
transitio n 20 ' 21 . Upon doping, the anomaly is suppressed 
monotonically and soon becomes an unremarkable kink 
at lower temperatures. No clear kink can be identified for 
x > 0.3. The suppression of this anomaly with P-doping 
is consistent with the suppression of AFM order of Fe 
ions, which is confirmed by a recent neutron study^. We 
mark arrows in Fig. 2(a) to show the Fe-AFM transition 
(T/vi, data from RefiSS). (iii) The resistivity upturn at 
low temperatures for x < 0.4. The upturn behavior is 
suppressed by doping or magnetic field, see the inset in 
Fig. 2(b). (iv) The metallic behavior in the entire mea- 
sured temperature range for x > 0.4. In particular, the 
resistivity drops rapidly at low temperatures for larger x. 
Such a drop in resistivity has already been observed in 
the HF metal CeFePO in the previous report^. However, 
no SC is observed down to T=2 K for the entire P doping 
range. 

It should be noted that the value of resistivity itself 
may not reflect the poly-crystalline sample's quality. Ac- 
tually the most likely impurity phases in the Ln-1111 
pnictides could be FeAs or Fe2As, and L112O3. Both 
FeAs and F82AS are highly metallic, more metallic than 
CeFeAsO itself at low temperatures. The metallic impu- 
rity phase like FeAs could distribute most likely in the 
grain boundaries and thus result in a lower resistivity if 
the main phase is a bad metal. The intrinsic resistiv- 
ity behavior of Ln-1111 parent pnictides below the Fe- 
AFM order is still ambiguous. For example, the resistiv- 
ity of single crystal of LaFeAsO shows an upturn at low 
temperatures^, meanwhile the resistivity of both poly- 
crystalline and single-crystalline CeFeAsO samples shows 
clear metallic behavior— 124. However, our more than 20 
CeFeAsi-^P^O samples with different P content were 
prepared under same conditions and the XRD patterns 
show that they all are of same high purity, the system- 
atic decrease in resistivity indeed suggests the increasing 
metallicity with increasing P content. 

The temperature dependence of dc magnetic suscepti- 
bility is shown in Fig. 3. In the high temperature range, 
x(T) increases with decreasing temperature down to 150 
K following a Curie- Weiss law for all < x < 1. The 
derived effective moment Heff is 2.54 /j,b for x = 0, and 
2.56 hb for x = 1 [insets of Fig. 3], very close to that 
of a free Ce 3+ ion, 2.54 fis- By lowering temperatures, 
a clear peak was observed at 4.16 K for x = 0, related 
to the formation of the Ce 3+ AFM order. With increas- 
ing x, x{T) increases and the peak becomes round (or 
shoulder- like) ; meanwhile the corresponding AFM order- 
ing temperature does not change too much. For x > 0.4, 
an obvious divergence between zero-field-cooling (ZFC) 
and field-cooling (FC) is seen at lower temperatures [see 
Fig. 3(b)], manifesting the FM ordered ground state of 
the Ce-4/ electrons. The suppression of AFM and devel- 
opment of FM are consistent with the enhanced metallic 
transport behavior in the course of P-doping. The FM 
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FIG. 1: (Color online) (a), X-ray powder diffraction (XRD) patterns of x=0, 0.2, 0.4, 0.6, 0.8 and 1.0 measured at room 
temperature, (b), Expanded XRD patterns showing the peak shift with the P doping, (c), Lattice parameters as functions 
of P doping level, (d-e), Rietveld refinement profiles of CeFeAsO and CeFePO. The refined atomic coordinates (x, y, z) are: 
Ce(0.25, 0.25, 0.1411), 0(0.75, 0.25, 0), Fe(0.75, 0.25, 0.5), As(0.25, 0.25, 0.6547) for CeFeAsO; and Ce(0.25, 0.25, 0.1508), 
0(0.75, 0.25, 0), Fe(0.75, 0.25, 0.5), P(0.25, 0.25, 0.6384) for CeFePO. 



ordered state persists until x ~ 0.9 where the divergency 
of x(^) vanishes. 

The evolution of Ce 3+ magnetism can be further 
demonstrated by isothermal field-dependent magnetiza- 
tion measurements, which are shown in Fig. 4. The tem- 
peratures are fixed at 2 K and 10 K, respectively. For 
x = and T=2 K, a kink in the magnetization curve 
is clearly seen, indicating a possible spin flop in the 
AFM state of Ce 3+ . When x increases the kink becomes 
weaker. As x > 0.4, M(H) starts to increase rapidly for 
small fields and tends to saturate at higher fields, i.e., 
it exhibits a nonlinear field dependence, see Fig. 4(d) for 
x = 0.4. For 0.5 < x < 0.9, a clear hysteresis loop can 
be observed and the saturated M value reaches the max- 
imum around x = 0.6 [see Fig. 4(c)]. Therefore, the long- 
range FM order exists in a wide doping range in between 
x ~ 0.4 and x <~ 0.9. The largest saturated magnetic mo- 
ment is about 0.95 (j,b, close to 1 us expected for a Ce 3+ 
doublet ground state^^. For 0.9 < x < 0.95 the hystere- 
sis loop is hardly observable down to the lowest measured 
temperature while the saturation tendency persists. This 
may be due to strong FM fluctuations at the quantum 



critical point associated with the 4 /-electrons. On the 
other hand, both the susceptibility and magnetization 
show a non-magnetic HF-like phase for x > 0.95 [see 
Fig. 3(b) and Fig.4(g-h)], consistent with the resistivity 
measurement. 

The metallic HF-like behavior can be further man- 
ifested by the specific heat measurement as shown in 
Fig. 5. For x = 0.95,0.97, and 1, the specific heat co- 
efficient 7(T) = C(T)/T shows logT behavior at low 
temperatures down to 2 K, see the inset of Fig. 5. We 
define a characteristic temperature T 7 (below which HF- 
like behavior appears) as the intersection of two extrap- 
olate lines of the specific heat data. For the x = 1 sam- 
ple, T 1 ~ 8.6 K. The saturated 7 should be about 700- 
800 mJ/molK or larger, in agreement with the result of 
Refi£ where 7 increases in proportion to logT before it 
gets saturated below 1 K. As a comparison, j(T) shows 
a A-peak maximized around 4.2 K and 6.8 K in the AFM 
and FM ordered phases, for x =0 and 0.8, respectively, 
in agreement with the magnetization measurement. A 
sudden increase in the saturated 7 close to 1 ~ 0.9 can 
be inferred by the values of j(T — 2K), accompanied by 
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FIG. 2: (Color online) Temperature dependence of resistivity 
in CeFeAsi-^P^O. (a), For < x < 0.4. Arrows denotes the 
temperature of Fe-AFM order, Tjvi, which were taken from 
Ref.— . (b), For 0.4 < x < 1. The inset shows the resistivity 
of £=0.3, 0.4 and 0.5, under zero field (solid lines) and (j,oH 
of 5 T (dashed lines). 



a decrease in the saturated magnetization under hqH = 
5 T. This implies the enhancement of the Kondo screen- 
ing near the FM instability. It should be noted that a 
very tiny kink was observed around 6 K on the C /T — T 
curve of x=0.95. We ascribe this additional peak to the 
contribution from some impurities with lower P content. 
Roughly estimated from the intensity of the kink, the 
amount of the impurity phase should be less than 2%. 

We summarize the experimental results by suggesting 
an electronic phase diagram shown in Fig. 6. At low tem- 
peratures, both the Fe-3d and Ce-4/ electrons show the 
long-range AFM ordering for x < 0.37. The d-electrons 
are Pauli paramagnetic(PM) for x > 0.4, while the /- 
electrons are FM ordered until x ~ 0.9. For x > 0.95, 
the /-moments are completely quenched and the whole 
system becomes a nonmagnetic HF-like metal. The HF- 
like region is roughly depicted by T 7 . Therefore, in ad- 
dition to the first magnetic QCP (denoted by x Cl ) asso- 
ciated with the AFM-PM transitions of the d- electrons, 
there could be a another QCP associated with the FM- 
HF transiton of /-electrons around x of 0.92 (denoted by 
x C2 ). x Cl has been known to be 0.4 by the study of neu- 
tron scattering^ 2 -. A turning point (denoted by x C3 ) as- 
sociated with the AFM-FM transition of the /-electrons 
at x C3 ~ 0.37 is quite close to x Cl . Meanwhile the strong 
FM fluctuations show up in the vicinity of x C2 . 




FIG. 3: (Color online) Temperature dependence of dc mag- 
netic susceptibility under ZFC (solid symbols) and FC (open 
symbols) protocols for CeFeAsi-^P^O samples under a mag- 
netic field of 1000 Oe. Arrows are guides to eyes to indicate 
the FM transition temperature Tc which is defined from the 
extrapolation of 1/x — T plot. The insets show the Curie- 
Weiss behavior of susceptibility above 150 K for CeFeAsO 
and CeFePO respectively. 



It should be noticed that the transition temperatures of 
the Fe-AFM ordering and the structure distortion can be 
approximately determined from the resistivity anomaly 
as shown in Fig. 2 and neutron studies^. They are undis- 
tinguishable in resistivity when x > 0.2. The neutron 
studj*22 revealed that both of them are suppressed at 
x Cl ~ 0.4, confirming an early theoretical prediction for 
the existence of a unique magnetic QCP driven by P- 
doping 25 . Moreover, the d-electron QCP seems to be very 
close to the turning point (x C3 ) of the /-electrons. This 
fact may indicate that the <i-electron magnetic transition 
is of weak first order accompanied by a reconstruction of 
the electronic states across x Cl . It is not certain if the 
coincidence implies the AFM ordering in Ce subsystem 
being supported by the AFM order in the Fe plane, be- 
cause the ordering patterns of the two subsystems arc 
different and the /-electron AFM order can coexist with 
the d-electron SC order in the F-doping case^^. Besides 
this puzzling feature, other implications of our experi- 
mental results should be addressed as follows. 

First, as the suppression of the a-axis is only a frac- 



5 




FIG. 4: (Color online) Isothermal magnetization of CeFeAsi-^P^O. (a), x = 0. (b), x = 0.2. (c), x = 0.3. (d), x = 0.4. (e), 
x — 0.6. (f), x = 0.8. (g), x = 0.9. (h), x = 1.0. For clarity, only the data of 2 K (solid lines) and 10 K (dashed lines) are 
shown. The insets of (e) and (f) display the magnetic hysteresis loop at 2 K. 
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FIG. 5: (Color online) Specific heat measurements of 
CeFeAsi-^P^O. A A-peak can been observed around T = 4.2 
K and 6.8 K for x — and 0.8, respectively, corresponding to 
the AFM and FM transitions. Inset: For x = 0.95, 0.97 and 
1, the specific heat coefficient j(T) = C(T)/T shows logT 
behavior at low temperatures down to 2 K, and tends to sat- 
urate to 7(0) when T — > 0. The 7(0) value was extrapolated 
to be about 700-800 mJ/K 2 mol or larger for x = 1.0. 



tion of that of the c-axis upon P-doping, the electronic 
states within the Ce-0 layer should not change signif- 
icantly. This is also supported by the non-magnetic 
LDA calculation from which the different Fermi sur- 
faces of CeFeAsO and CeFePO can be mainly ascribed 
to the pnictogen heigh t 18 ' 26 . Therefore the HF behav- 
ior in the end compound, CeFePO, is originated unam- 
biguously from the inter-layer coupling between the Fe- 
3d and Ce-4/ electrons. Second, the /-electron AFM- 



FM transition could be explained by assuming the ex- 
istence of a bare d-f Kondo coupling (Jk) in the or- 
dered phases, as it can mediate the indirect Rudcrman- 
Kittcl-Kasuya-Yosida (RKKY) interaction (Jrkky) be- 
tween the Ce-4/ local moments. Intuitively, the exchange 
interaction between the local /-spins , Jf, is given by 



Jf = J^p + Jrkky, with jy 1 being the exchange in- 
teraction in the absence of Jk, and | Jrkky\ ~ Jr^p, 
with Np being the density of states at Fermi level (of 
mainly <i-clcctron characteristic). Jrkky is negative if 
Np is small, as it should be since the parent compound 
CeFeAsO is a bad meta l 12 ' 13 . Upon P doping, both Jk 
and Np should increase, leading to a sign change in Jf. 
Third, the negative Jf is also the driving force for the 
strong FM fluctuations in the vicinity of the critical point 
x C2 , where the Kondo scale Tr (which increases exponen- 
tially with JrNr) starts to dominate over the magnetic 
ordering energy scale, Tf ps \Jf\. 

It is well-established that the RKKY interaction be- 
tween the local moments is mediated by the Kondo- 
coupling between the itinerant charge carries and local 
moment a 27 ' 28 . The Kondo coupling can lead to both the 
Kondo effect and the RKKY interaction, with the latter 
being responsible to the magnetic ordering of the local 
moments. The crucial point is that the dependences on 
the Kondo coupling (together with the DOS of the charge 
carries) of the Kondo energy scale (describing the Kondo 
effect) and of the RKKY interaction (characterizing the 
ordering temperature of the local moments) are differ- 
ent, as described by the intuitive formulas in the inset 
of Fig. 6. Based on the similar intuitive formulas, Do- 
niach first suggested that a quantum phase transition 
from magnetic ordering phase to the paramagnetic HF 
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FIG. 6: (Color online) Electronic phase diagram of 
CeFeAsi-^P^O (0 < x < 1)). The shaded region (in red) 
in the left side shows the antiferromagnetic ordering of Fe- 
3d moments (denoted by Fe-AFM below Tjvi— ). In the lower 
temperature portion of this region, simultaneously, Ce-4/ mo- 
ments order antiferromagnetically (marked with Ce-AFM) 
below Tat2- The small pink area at the lower right cor- 
ner displays the nonmagnetic state with heavy fermion (HF) 
behaviour. In the lower middle regime (in green), the 4/- 
moments are ferromagnetically ordered (labeled by Ce-FM) 
below Tc- Inset: The schematic competition between Jj and 
T K accounting for the the AFM-FM-HF evolution of the /- 
electrons. 



AFM ordering gap of the d- electrons^. This picture pro- 
vides a natural explanation for the strong FM fluctua- 
tions in a class of homologous HF compounds including 
CeFePO^ and CeRuPC^ 30 , Therefore the observed f- 
electron phase diagram of CeFeAsi-^P^O system is con- 
sistent with the Doniach picture^!. In addition, it also 
predicts a possible weak Kondo phase in a narrow doping 
region enclosing the turning point x Ca (at very low tem- 
peratures). As x Ca is very close to x Cl , it may account for 
the observed suppression of the SC near the AFM QCP. 
The existence of this peculiar phase and its possible con- 
nection with the transport and magnetic properties still 
need to be clarified. 

IV. CONCLUSION 

In summary, our experimental study on 
CeFeAsi-zP^O (0 < x < 1) reveals a rich phase 
diagram which consists of an AFM quantum phase 
transition of the Fc-3rf electrons and an AFM-FM-HF 
evolution of the Cc-4/ electrons. The AFM QCP of 
the d-electrons is very close to the turning point of the 
AFM-FM transition of the /-electrons, while the d-f 
hybridized HF state is separated by an /-electron FM 
instability. In contrast to the F-doping case where the 
high-T c SC is induced, no SC is observed down to 2 K in 
the present case, signalling the unique role of P-doping 
in suppressing the d-electron correlatio n 13 ' 25 . These 
results highlight the importance of the inter-layer Kondo 
physic a 15 ' 31 and the interplay of the d- and /-electron 
correlations in the rare-earth iron pnictides. 



phase could exist in a Kondo lattice^!. Here, we sug- 
gest that the AFM-FM-HF evolution of the /-electrons 
in CeFeAsi-zP^O is due to the competition between Tk 
and Jf (instead of Jrkky in the conventional Kondo 

lattice^), see the inset in Fig. 6. Notice that J^p is 
positive and relatively small in the absence of Jrkky- 
Thus in the AFM phase Tk is mainly suppressed by the 
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